Effective Hamiltonian and low-lying energy clustering patterns of four-sublattice 
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We study the low-lying energy clustering patterns of quantum antiferromagnets with p sublattices 
(in particular p = A). We treat each sublattice as a large spin, and using second-order degenerate 
perturbation theory, we derive the effective (biquadratic) Hamiltonian coupling the p large spins. In 
order to compare with exact diagonalizations, the Hamiltonian is explicitly written for a finite-size 
lattice, and it contains information on energies of excited states as well as the ground state. The 
result is applied to the face-centered-cubic Type I antiferromagnet of spin 1/2, including second- 
neighbor interactions. A 32-site system is exactly diagonalized, and the energy spectrum of the 
low-lying singlets follows the analytically predicted clustering pattern. 
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I. INTRODUCTION 

Many frustrated isotropic antiferromagnets develop 
long-range order, but their classical ground states have 
nontrivial degeneracies, so that the true (quantum) spin 
ordering pattern is not obvious and may be decided by 
small perturbations. Exact diagonalization calculations 
have been carried out on finite systems and the low-lying 
eigenenergies show interesting clustering patterns. In this 
paper we develop an effective Hamiltonian approach to 
studying energy clustering for a finite-size lattice. The 
approximation that we use is to divide the system into 
sublattices following the classical ground states, represent 
the spins in each sublattice by one large spin, and ac- 
count for the fluctuations within sublattices by an effec- 
tive Hamiltonian, which couples the large sublattice spins 
to each other. This effective Hamiltonian is much more 
manageable (as far as exact diagonalization is concerned) 
than the original Hamiltonian. From this effective Hamil- 
tonian, using notions of Bohr-Sommerfeld quantization, 
tunneling, and the geometrical phase, one can predict 
the energy patterns of low-lying singlet states of a small 
system, as may be "observed" in exact diagonalizations. 
Recently, in Ref. |], two of us studied the triangular anti- 
ferromagnet with nearest and next nearest neighbor in- 
teractions and compared analytical predictions with ex- 
act diagonalization results in Ref. ^. In this paper, we 
present a systematic derivation of the (biquadratic) ef- 
fective Hamiltonian for a class of antiferromagnets (this 
was explained in Ref. |l| merely by a footnote). And we 
check its predictions on the eigenenergy clustering pat- 
terns by carrying out exact diagonalizations of the fee 
Type I antiferromagnet. We will also correct some nu- 
merical mistakes in a table in Ref. |l|, resulting in a much 
better agreement between diagonalization and analytical 
results for the triangular lattice. 

We start with a general Heisenberg Hamiltonian 



where i and j run over N lattice sites and each has 
quantum spin s. Jij is the coupling constant between 
the sites i and j and is equal to J„i when i and j are 
rn-th nearest neighbors (m = 1, 2, ...). 

We assume that the ground state of H, has long-range 
antiferromagnetic order (of the spin directions). In the 
simplest antiferromagnets (such as the square or triangu- 
lar lattice with nearest-neighbor interactions), the clas- 
sical ground state is unique (modulo rotations) and the 
quantum ground state has the same ordering pattern. 
The antiferromagnets of interest here have continuous de- 
generacies - not due to symmetof_- which are typically 
broken by quantum fluctuationsjj'u 

The present work applies to a subset of all classically 
degenerate Heisenberg antiferromagnets, which we call 
"p-sublattice." This means the lattice divides into p 
equivalent sublattices, and the classical ground states are 
precisely those in which (i) the spins are parallel within 
each sublattice; (ii) the vector sum of all the sublattice 
magnetizations is zero. The two main examples that be- 
long to this class and have motivated this study are the 
triangular and fee Type-I systems with appropriate Ji 
and J2 interactions. The p-sublattice antiferromagnets 
are discussed in more detail below (Sec. ||). 

The effective Hamiltonian Tdcs is a function of the p 
sublattice spins. The simplest terms of the required sym- 
metry are of biquadratic form, and any more complicated 
form would be difficult to use in an analytic calcula- 
tion. This effective Hamiltonian is derived by expanding 
around mean-field theory. That is, one invents a "mean- 
field Hamiltonian" TYmf which has the same classical 
ground state as the true one, and which has an exactly 
known quantum ground state. The full Hamiltonian is 
then written as 



(1.2) 



(1.1) 



and TLcS is obtained by second-order perturbation theory 
in Sii. 

Mean-field theory may be set up in two different ways 
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and these produce two different recipes for Ti-cSn The 
first approach was foUowed W Larson and Henleju, gen- 
erahzing a formula of Longcl. There, Hmf consists of 
local fields that fix each spin in its classical direction; the 
result is a derivation of Tioft, usable for any antiferromag- 
nct with a nontrivial degeneracies, but only for an infinite 
system. In this paper, we follow the second approach to 
mean-field theory: that is, the zeroth-order Hamiltonian 
is an infinite-range model in which the individual spins 
and their products are replaced by the sublattice aver- 
ages. This has the advantage that the finite-A^ effects 
are accounted for (this is important for comparing with 
exact-diagonalization results). 

The calculation proceeds in several stages. The first 
stage is to map each sublattice in the iV-spin systeip ap- 
proximately to just one spin of maximum lengthB The 
net energy from the neglected spin fluctuations of the 
spins in each sublattice (away from the perfectly aligned 
state), is approximated by an effisctive biquadratic inter- 
action favoring coUinear states E'LJ Such interactions ap- 
peared first in the theory of selection among degenerate 
classic^rfflound states (the so-called "order-by-disorder" 
effect).™ 

Mathematically, our calculation involves two steps. 
First, operator equivalents that operate on sublattice 
spin states are obtained for one spin and the product of 
two spins. This step is purely mathematical and follows 
from the Wigner-Eckart theorem concerning the addition 
of angular momenta. The second step is summing the op- 
erator equivalents over the lattice, and this step depends 
on the geometric relationship of the antiferromagnet's 
sublattices. 

In this paper, we derive the biquadratic effective 
Hamiltonian and apply it to the fee Type I antiferro- 
magnet (with nearest- and next-nearest-neighbor inter- 
actions) . We generalize this derivation to treat arbitrary 
(translation-invariant) exchange couplings, in particular 
to properly handle couplings which connect two spins of 
the same sublattice (such as J2 in the fee Type I anti- 
ferromagnet) . Combining this with the remaining results 
of Ref. |l| on Bohr-Sommcrfeld semiclassical quantization 
and geometrical phase, we find that the expected split- 
ting pattern for the fee Type I indeed agrees with numer- 
ical diagonalization results for TV = 32 sites. (See Ref. ^ 
for a more detailed description of numerical diagonaliza- 
tion of fee systems and Ref. |l^ for an earlier study of 
another spin system, a three-spin cluster, using many of 
the key semiclassical ideas here.) 

This paper is organized as follows: we first introduce 
the class of p-sublattice antiferromagnets and give some 
examples (Sec. ||); we then motivate and give the mean- 
field Hamitonian (Sec. [H); next, the biquadratic effec- 
tive Hamiltonian is derived from second-order pertur- 
bation theory (this is the core of the paper) (Sec. IV). 
We then apply semiclassical considerations to the result- 
ing effective Hamiltonian (Sec. 0) analytically calculat- 
ing the energy-level pattern. The analytic result for the 
type-I fee antiferromagnet, as well as a convenient four- 



spin exact diagonalization based on the same effective 
Hamiltonian, is compared with exact diagonalization of 
a 32-site system ( Sec. JVj ); finally, we make some closing 
observations (Sec. [VI I| ). 



II. THE P-SUBLATTICE ANTIFERROMAGNETS 

We will set up a formalism that includes all realizations 
of p-sublattice antiferromagnets. It is convenient to de- 
fine a geometric factor z™^, as the number of m-th neigh- 
bors on sublattice a of a site on sublattice (3. (We will 
consistently use Greek indices a,P to label sublattices 
and /i, v, p to label spin components.) Notice that is 



not an attribute of the Hamiltonian (1.1), but of the ge- 
ometrical arrangement of sublattices within the system. 

Unless specifically noted, we will assume (and need) 
sublattice-pair permutation symmetry. This means that 
if we perform any of the p\ permutations of the sublattice 
labels, the geometric relationship between any pair (a, (3) 
will be unchanged (modulo lattice symmetries such are 
translations, rotations, and reflections). 

The classical ground state manifold is labeled by the 
unit vectors {rhc}, a — 1, ... ,p for the directions of the 
sublattices. Since we seek a quantum-mechanical Hamil- 
tonian, however, it should be expressed in terms of the 
corresponding operators, which (apart from length nor- 
malization) are the sublattice spins. 



(2.1) 



i:a(i) — a 



where a{i) — 1,2, gives the sublattice index of the 
spin Si. Here the sum is over N = N/p sites within the 
sublattice a. Each can be in states with total spin 
ranging from zero to = Ns. The classical state corre- 
sponds to a coherent state in which every in sublattice 
a is aligned with riia, hence each sublattice spin has its 
maximum length S. 

Our final goal will be an effective Hamiltonian Tioff de- 
fined on the sub-Hilbert-space of "p-sublattice states", 
being all states in which every sublattice spin has maxi- 
mum length, = S{S+1)] thus Tioff will have the form 
of a spin Hamiltonian. (Note that these "p-sublattice 
states" have maximum length. The sublattice magneti- 
zation is not determined and will be decided by diagonal- 
izating the secular Hamiltonian matrix as developed in 
the degenerate perturbation below.) Along the way, we 
will often deal with the Hilbert space of one particular 
sublattice - e.g. the operator equivalents are computed in 
this space. In particular, we let \(j)a) mean any eigenstate 
with — S{S + 1) (maximum length) and simultane- 
ously a definite S^. The direct products of such states 
make up the basis for p-sublattice states. On the other 
hand, \tpa} will be used for a simultaneous eigenstate of 
and with any total spin, < Sa < S. 

The chief concrete examples havfl-H — 4, being: (i) the 
J1-J2 triangular antiferromagnet with Ji > and 
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second-neighbor coupling J2 € (Ji/8, Ji) - the effective 
Hamiltonian for this case was latmpunced in Ref. Q (ii) 
the Type I fee antiferromagnetll3ll3 with Ji > 0, J2 < 0. 
In either example, there are only nearest and next- 
nearest neighbor interactions (J„i = for m > 2); the 
Type I fee differs from the triangular case in that the 
next-nearest-neighbor couplings (z^) connects sites of the 
same sublattice. (Table |, below, includes the values 
for these lattices.) The honeycomb and pyrochlore lat- 
tices also support four-sublattice states for appropriately 
chosen exchange couplings (Ji and J2). After allowing 
for rotations, a 4-sublattice antiferromagnet stili iias a 
two-parameter family of classical ground states.ll2nl3 

The case p = 3 is less interesting: there is just one 
way to add up three equal-length spins to make total 
spin zero, so the classical or "mean-field" ground state 
is already unique. But our approach would be applica- 
ble to the three-sublattice antiferromagnet in an external 
field, in which case the ground state does have nontrivial 
degeneracy. 

There exist p-sublattice antiferromagnets with p > 4, 
and possessing "sublattice-pair permutation symmetry" , 
although their Hamiltonians are physically improbable. 
For example, the p = 5 case is realized by the x ^/E 
sublattice; in a square lattice without reflection symme- 
try. The ground state has five sublattices, e.g., when 
the couphngs to the ±[2,1] and ±[1,— 2] fourth neigh- 
bors are strongly ferromagnetic, and the nearest neighbor 
couplings are antiferromagnetic. A similar construction 
with p — 7 may be made using the ^/^ x sublattice of 
a triangular lattice. 



III. MEAN-FIELD HAMILTONIAN 

T he o bjective of this section is to construct from H of 
Eq. (1.1) a "mean-field" Hamiltonian T^mf that is defined 
on the entire Hilbert space but expressible in terms of the 
sublattice spins {Sq,}. In other words, it is a projection 
of the exact Hamiltonian onto a kind of infinite-range 
Hamiltonian. We will show that the ground states of 
Hmf are all the p-sublattice states; this high degeneracy 
corresponds to the classical degeneracy. The real val ue o f 
TYmf is that we can split the Hamiltonian in form (1.2) 
and obtain Tieff via second-order perturbation in STi. (see 
Sec. 0). 

The eigenstates of the lattice system that correspond 
one-to-one to those of our p-spin effective Hamilto- 
nian are not literally superpositions of the "p-sublattice 
states." The true eigenstates include quantum fiuctu- 
ations that admix other states (as is implicit in our 
second-order perturbation), so the overlap with the p- 
sublattice subspace must decrease exponentially with in- 
creasing system size N; indeed the autiiDrs of Ref. |2| com- 
puted the overlap and found it smallJlil The low-energy 
states of the lattice antiferromagnet are related to the 
p-spin system's eigenstates much as the excited states of 



a Fermi liquid are related to those of noninteracting elec- 
trons.Ej In this paper, our main interest is in analyzing 
the clustering patterns of the low-lying energy spectrum. 
We will focus on calculating eigenenergies, not magneti- 
zation expectations. 

The mean-field Hamiltonian TCmf, to be a function 
only of the sublattice totals {Sq}, must be symmetric 
under any way of permuting the spins within a sublat- 
tice. Furthermore, the states that are symmetric under 
such permutations are just the p-sublattice states, and 
Hmf ought to have the same expectation as H in such a 
state. (This will be demonstrated in the next section.) 

The unique way to project 7i onto Hmf and satisfy 
these conditions is to average each term Sj • Sj over per- 
mutations of the sublattice(s) to which the sites i and 
j belong. There are two cases. First, the two spins are 
on different sublattices {a(i) 7^ ct(j)); then we simply re- 
place s^-Sj {Sa(i)/N)-{Sa(j)/N). Second, the spins are 
on the same sublattice (a(z) = a{j) = a with i ^ j); in 
this case there are N{N — 1) pairs of such scalar products 
within one sublattice. (No permutation can ever make 
i — j.) We use the identity = Ns{s -I- 1) -I- Si/j ' 
to obtain that • Sj (S^ -iVs(s + l))/iV(7V- 1) in this 
case. (When Sa = S, i.e., maximum length, this average 
is simply s^; restricted to the maximum-spin subspace, 
the projection we are doing here agrees with the opera- 



tor equivalent of • Sj, defined in (4.6), below.) Finally, 



putting these two kinds of terms together and performing 
the sum, we get 



Hmf — 2 ^ Jap^a ■ 



Ns(s + l' 



where 



^ ^(J™Z™^)/7V. 



(3.1) 



(3.2) 



We restrict our attention to the class of p-sublattice 
antiferromagnets with sublattice-pair permutation sym- 
metry, so = for each m and a ^ (3 and z™^ — . 

Thus = J12 for slU a ^ l3 and Jaa = Jii for all a's. 
The mean-field Hamiltonian becomes. 



Hmf - 2^12 (^lot - X! 



(3.3) 



1 1 1 

= 2^i2SL ~ E Si - 2P'^ii^^(^ + 1)' (3-4) 

a 

where Stot ~ J2a and we have also defined 
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Js — Ji2 — 



N 



N -1 



Jii 



(3.5) 



Next we consider the ground states for Hmf, which we 
denote |<i>). For this, we further restrict to the situation 
wher e Js > 0, i.e., the coefficient in front of ^ in 
( |3.4D is positive. (This condition is satisfied by both the 
triangular antiferromagnet (Jn = 0) and the fee type-I 
case ( Ji2 > and Jn < 0). See Table | for a collection 
of relevant parameters.) 

Then |<i>) should obviously satisfy StoJ*) = and 
S^|$) ^ S{S +1)1$) for every a. In other words, the 
ground states of Hmf are just the p-sublattice states 
which are also total singlets, so we call them "p-sublattice 
singlet states." (A similar result is obtained for the clas- 
sical ground state: that each of the p sublattice spin vec- 
tors is at maximum length and their vector sum is zero.) 
For p — 4:, there are 2S+1 such singlets;! for p = 5 and S 
is an integer, there |are 5S{S+l)/2 + l (see Appendix^). 
Lecheminant et aln already recognized that, in the Ji- 
J2 triangular antiferromagnet, the lowest energy singlets 
correspond to the four-sublattice singlet states. 

We emphasize here that we have a degenerate family 
of singlets. These states have the rotational symmetry as 
the true finite-N ground state is known to have (i.e., they 
are singlets). The sublattice spins have maximum length 
but the z components (magnetization expectations) are 
not determined. Which linear combination of these states 
best represents the ground state is only settled by the 
degenerate perturbation calculation that follows, which 
diagonalizes an effective Hamiltonian within this degen- 
erate family of singlet states. 



The ground state energy of (3.4) is 



Emf = -^p 



(3.6) 



Table | gives Emf for the triangular and fee lattices. 

It is necessary to correct this mean-field Hamiltonian 
by ficS, the effective biquadratic Hamiltonian,EI to pre- 
dict the sequence of quantum numbers of the low-lying 
states, and the clustering pattern of their eigenenergies. 



IV. BIQUADRATIC EFFECTIVE HAMILTONIAN 

The calculation of the second-order effective Hamil- 
tonian TiesiiSa}) will be presented now, and it is the 
central analytic result of this paper. We will use the 
second-order degenerate perturbation theorjtj for ST-i = 
H — Hmf, 

(.'i„,.i.).-^(£M|i>«, ,4,1, 

where 1$) and |<I>') are both p-sublattice singlet states 
(i.e., in the ground state manifold of TYmf) and j^*) is an 



excited state of T^mf- In the process, we will also show 
that this choice of mean field Hamiltonian enables the 
first-order correction to vanish, i.e.. 



It is easy to show that 
where, 



(4.2) 
(4.3) 



m-th(ij) 



^E^™ E s,.s,-i^z™,S„.sJ , (4.4) 



-th(i3) 



and, 



-th(ij) 

= a(j) 



(4.5) 



and Ssi — Si ~ Sa{i)/N. Here the subscript "diff" means 
that the sum is over {ij) pairs on different sublattices 
and "same" means that they are on the same sublattices; 
also, "m— th(ij)" means a sum is restricted to pairs 
which are m-th neighbors. We should consider what kind 
of excited states l^") are present in the sum of the second- 
order perturbation theory (4.1). Because this consider- 
ation in part relies on the Wigner-Eckart theorem and 
the concept operator equivalent, we will turn to these 
mathematical concepts now. 



A. Operator equivalents 

The general problem of this sort of calculation is to 
take an operator defined in a large Hilbert space, restrict 
this Hilbert space to a smaller one, and re-express its 
action by an operator acting on the restricted space. The 
latter operator will be called the "operator equivalent" 
of the first one. In the present paper, we start from an 
operator 0{{si}) that acts on all states of the spins {si} 
belonging to sublattice a and convert it into an operatpt 
6(Sq) acting on states of a single net sublattice spin Sq.EZI 

The obvious [-tool for this calculation is the Wigner- 
Eckart theoremEEl which deals with angular momenta ad- 
dition. It says that among states of a single spin multiplet 
(such as the maximum-spin states |0q)), any operator 
that transforms as a particular representation of the rota- 
tion group, has matrix elements given by Clebsch-Gordan 
coefficients, apart from a single constant factor setting 
the overall scale. Hence we can replace that operator by 
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a simpler one belonging to the same representation: here, 
a polynomial 6(Sq,). There is a unique such polynomial 
for each of the rotational symmetries (scalar, vector, or 
traceless-symmetric tensor) . Specifically, we will be given 
operators 0{{si}), with a given symmetry, and write 
their matrix elements {(j)'^\0{{si})\(j)a) = c{(j)'^\d{Sa)\4>a) 
where c is a constant that depends on Sa, and possibly 
on i, but not on the z-component of the spin (see e.g., 
Ref. pj| ). In our calculation, we restrict our attention 
to the maximum spin manifold Sa = S because it turns 
out that we will need only matrix elements between the 
quantum p-sublattice singlet states |<i>) and |$'), where 
Sa = S for all a; in any case, when Sa < S there is more 
than one multiplct and the Wigner-Eckart theorem is not 
enough to specify the matrix elements. 

The results for scalar and vector operators are. 



'l0o)=O, 



(4.8) 



a{i) = a{j) = a, i 7^ j, Sa^ S (4 



1 



a{i) — a, Sa — S. 



(4.7) 



where the l ast s tep relies on the scalar operator equiva- 
lent result (16). This tells that SHsamc gives no con- 
tribution to the second-order perturbation (O), i.e., 
(^|(57isamc|*i') = 0. Of course, the argument here also 
shows that (57isamo gives no contribution to the first-order 
perturbation (|4^), i.e., (^'I^Tisamd*) = 0. 

In the same way, we turn to the contribution by (57idiff • 
It is easy to show that [{Ssi ■ Ssj), S^^J = which gives 
[(5?i:diff,S?ot] = 0- Since |$) has S^^^t = 0, it follows that 
S^ot = for I*) also, otherwise (^'I^Hdiffl*) would be 
zero. 

Furthermore, that matrix element can be broken 
into terms of form {^piip2''p3'>p4\^^i ' '5sj|0i020304). As- 
sume without loss of generality that a{i) — l,a{j) — 
2, so the matrix element factors as {ijji\dsi\(l)i) ■ 
(?/'2|i5sj|02)(^3|03)('04|04)- Consider the first factor: as 
6) defined above, \4>i) has Si = S, and the raising and low- 
ering operators in Ssi can change Sa by -1-1, 0, or —1 only. 
An increase 5 ^ S" -I- 1 is obviously impossible since S is 
the maximum spin. Also, when has the same length 



spin, S, the matrix element is zero by the construction 



The constant coefiicients, c = for the scalar operator -^f ' = S<,(,)/7V, as found using the operator equiv- 



equivalent and c = 1/N for the vector operator equiva- 
lent, were found by explicitly calculating the matrix ele- 
ments on both sides for a particular choice of states: \(j)a) 
and \(f)'^) are both the state with all spins aligned along 
+z, i.e., Sa = S. This is the most convenient choice to 
simplify the algebra. 

As noted above, every maximum-spin state \(j)a) is 
symmetric under permutations of the sites of sublattice 



hence our operator equivalent expressions (4.6) and 



alent (4.7), Hence, by elimination, the only states |-0i) 
that can give any nonzero matrix element in (ipilSsilc/ji) 
must have spin S — 1; the same holds for \'ip2), whereas 
iV'a) and IV'4) obviously have the same spin S" as do | (/is) 
and 1^3). 

Finally, every term of SHdiS contains products of two 
such operators 6si ■ Ssj, so {'^\6'Hdis\^} is nonzero only if 
1 5*) has exactly two of the p spins reduced by one. This 
gives 



(4.7) cannot depend on i and j. This symmetry could 
be used in pla ce o f the Wigner-Eckart theorem to derive 
eqs. ( Li ) and ( i/7), but only the Wigner -Ecka rt approach 
works in the tensor operator case (Sec. [VD, below). 



= -^JsiiP - 2)SiS +1) + 2(3 - 1)5) 
1 ~ 



(4.9) 



B. Energy denominators 

Now we are ready to compute the matrix elements 
(^'|(57tI|<I>) and denominators E^i ~ Emf for the excited 
states in the perturbation sum (4.1). First we address 
SHsamc which can be considered as acting on states of 
sublattice a, \(j>a) or \ipa)- Notice that the operator Sj__Sj 
with a{i) — ce{j) = a commutes with S^. From (4.5), 



therefore, [SHsamc, S^] = for all a. This means that, if 
the states [(pa) and {ipa) have different spin, the matrix 
element {ipa\SHsamc\4'a) is zero. Since in our perturba- 
tion sum 1$) is a maximum spin state, the state \(j)a) has 
spin S and so S'^ = S* for {ipa) also. Finally we have, for 
a{i) = a{j) = a, 



Subtracting the ground state energy Emf (|3.6| ), the en- 
ergy denominator in the perturbation sum ( [1.1| ) that can 
give nonzero contribution is always 

E^ - Emf = 2JsS (4.10) 

This argument also gives that \5Hdis\^) = within 
the ground states of T^mf- Together with the result 
(<i>'|(57isamc|'I') = obtained earher in this section, we 
conclude that the first-order correction to Hmf by SH is 
zero, i.e., eq. (4.2) is established. 



Sl-Ns{s + 1) 
N{N- 1) 



C. Reduction of Perturbation Sum 

We have now concluded that all the states j'l') that 
cont ribute nonzero numerator to the perturbation sum 
(4.1) have the same denominator. This key fact enables 
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us to factor the denominators out of sum (4.1). Then 
we can extend the sum to run over all states j^*), not 
just the excited ones. (Every added term is zero, since 
1 5") is a maximum-spin state in the added terms, and we 
just showe d th at (^'|(57i|$) is always zero in that case.) 
Equation (4.1) then becomes 



($'|7T!eff|$) = - 



2JsS 



(4.11) 



S o in stead of trying to sum over excited states in 
eq. (4.1), our goal now will be to find the operator equiv- 
alent of {SHdis)^ as in eq. (4.11). Squaring (4.4) gives 



{<^'\{6Hdisf\<f) = J™J„0™„|$), (4.12) 

mn 

where, 



A is a scalar operator, Vp a vector operator, and T^i, a 
second-rank symmetric traceless tensor operator. 
By the Wigner-Eckart theorem as stated in Sec. |IVj 



the operator equivalents of the three operators in ( 4.16| ) 
are known functions of Sq, apart from coefRcients de- 
pending only on (i,fc). As before, we obtain the coeffi- 
cients using the most convenient choice of states = S 
(all spins aligned up) to simplify the algebra. With 
a{i) = a{k) ~ a, for the scalar operator A, we obtain 



(4.17) 



Note we can have i — k here, because i and k come 
from two separate scalar products (s^ •Sj )(sfe -s;), whereas 
earlier in (^^) we necessarily had i ^ j. Here Sik appears 
because when i = k we have = s(s -I- 1). Following the 
same procedure for the vector and tensor operators (using 
all-spins- aligned states to compute coefficients), we get. 



m-th(ij) Ti-th(fcO 



{Si ■ Sj){Sk ■ Si) 



^ ^ S„ • S;3 • 



7V2 



(4.13) 



In deriving this, we have used another operator equiva- 
lent relation 

KKS:\K) = l{^'JS^S:;,\(bc.), a{z)=a, (4.14) 

which is true because being a state of maximum 

Sa, has the symmetry of permutation among all sites i 
in sublattice a. 



D. Operator equivalent of quadratic operator 



To obtain TYmFj we have considered, in section III, the 



operator equivalents of scalar and vector operators. Here 
for HcS , we need the operator equivalent of (s^ • s j ) (s^ • s; ) 
which can have two of the spin operators (e.g. and s^) 
on the same sublattice. For this case we will need the 
operator equivalents for tensor operators. 

As before, to find the operator equivalent we consider 
a state of just one sublattice |0q) with the maximum spin 
S. Recall the decomposition of a direct product of two 
vector operators: 



(4.15) 



2N 
1 2s 



iV 25*- 1 



(4.18) 



I^Z-a). (4.19) 



-(5^^5^ + 5^5^i)--V^(5 + l) 

And putting the three terms together as in ( [4.15| ), we 

get. 



1 2s - S^k 
N 2S-1 



+ Siks)Spt, 

Ja'-'a) 



— a{k) — 



(4.20) 



It may be c hecke d that summing both sides of ( |420| ) 
indeed gives (4.14). As another check, notice that when i 
and k are distinct sites belonging to the same sublattice 
a, the spin operators sf and s^ commute. Hence the 
operator equivalent of s^'sj^ must be symmetric under 
exchange of « ^ k, even though [S^,S^] ^ 0; indeed 
( 4.20D has this property. 

We return to consider matrix elements between p- 
sublattice states such as |$) and 1$'). To compute the 
effective Hamiltonian due to (STCdis)^, we w ill ne ed the 
operator equivalent of (s^ • Sj){sk ■ sj) (see ( 4.13D ). We 
give here the representatives of the three possible situa- 
tions. (All others are equivalent, by permutation of the 
sublattice indices.) 



1. Case a{i) = 1, a{j) = 2, a(fc) = 3, a{l) = 4 



where 



A — Si ■ Sfc, Vp — -Spiyp SiSk, 



T,,^-{stsl + st.4)--5p,A. 



(4.16) 



b'\{s, ■ s,)(sk ■ siM = (<?i'|^(Si • S2)(S3 • 



(4.21) 
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2. Case a{i) = a{k) = 1, a{j) = 2, a{l) = 3 



2. Case Q(i) = a{k) = 1, a{j) = 2, q(0 = 3 



y|(s,.s,)(sfc-sO|0) = (</.'||4^(S2-S3) 



G(Si, S2, S3) + -^(82 • S3) - G(Si, S2, S3) 



where 



ljSi.(S2XS3)](5.fc}|(/.), 



(4.22) 



G(Si, S2, S3 



1 



iV3(2S'- 1) 



-((Si.S2)(Si.S3) 



-(Si-S3)(Si.S2))-i^(^+l)S2-S3 



(4.23) 



Nz^Nz^ term s and Nz]^z^ of them have i = k (for 
the Sik terms in ( 4.22| )). The contribution to the total 
sum turns out to be simple 



-^zr2^r3(Sl-S2)(Si.S3). 



(4.27) 



We can obtain this resul t either from the complicated 
operator equivalent ( 4.22| ) or from the fact that here the 
double sums X)m-th En-th summed separately 

to give the sublattice spins. 



3. Case a{i) = a{k) = 1, a(j) = a{l) = 2 



3. Case a{i) = a{k) = 1, a{j) = a{l) = 2 



y|(s, • s,)(sfe • siM = + i^(Si, S2) 

is3-2sF(Si,S2)]((5,fc + <5,7) 



+ F(Si, S2) - ^^^(Si • S2)] (<5.fe^,7)} 



where, 



F(Si,S2) - ^ 



2Ar2 
1 



(4.24) 



(Si-S2)2 



iV2(25'- 1)2 

\si.S2)-i5^(^+l)^' 



(4.25) 



2' " 3 
These three formulas are obtai ned from summing over 
Ai, 1^ in E(sfs^)(SfcSr) and using (|3) and (^ ) for the 
operator equivalents of one single spin and the product 
of two spins on the same sublattice. 



This is the most complicated case. There are 
Nz'^Nz]^ terms and Nz]^z]^ of them have i = k (for 



the S,k term in ( 4.24[ )), Nz^z^ of them have j = I (for 
dik term), and Nz"2 of them have i = k and j = I (for 
Sik^ji term). Here the contribution to the total sum is 



^izT,ns,-s,f 



(Si • S2) 



(25-1)2 

5^(5-2) 
(25-1)2 



(Si.S2)2-i ( 1-^^ 

^ ^ 2 V (25-1)2 



+5mnNzT2 ( + ^F(Si, S2) - ^(Si • S2) ) , (4.28) 



where F is given above in (4.25). 



F. Final form of effective Hamiltonian 



E. Summing over the lattice 

With the operator equivalent of (s^ • Sj){sk ■ s;) we can 
sum over the lattice and obtain the operator equivalent of 
Oynn in ( |4.13| ). The results depend on how the sublattices 
are related geometrically to each other, via defined 
in Sec. ^ Again the representatives of the three possible 
situations are the following. 



1. Case a[i) = 1, a{j) = 2, a[k) = 3, a{l) = 4 

Here we have Nz^ pair of {a{i) = l,a(j) = 2) and 
TVzJJ pair of (a(fc) = 3, a{l) — 4) i n the double sum 
Sm-thSn-th- Contribution of ( 4.21 ) to the total 

sum is 



i-zr2^3"4(Sl-S2)(S3-S4). 



(4.26) 



Combining the three terms ( [4.26y4.27] , [4.28|) and their 
cyclic permutations, we get 

((S„ • Spf - 25(5 - 1)(S„ • S^) + 5^(5 - 2))! \^) (4.29) 



Finally, we perform X^nm ( 4.12| ) and use ( 4.11 ) and 
the fact that a singlet state |$) is simply a linear combi- 
nation of p-sublattice states |(/)) to obtain our result for 
the biquadratic effective Hamiltonian, 



Tic 



where 



•s, 



4JsiV5(25-l)2^^^l 
25(5 -1)(S„-S;5) + 5^(5 -2) 



(4.30) 
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K^p^J2(-J,-)''»0-j(j2jmZ^p] ' (4-31) 

m \ ni / 

and we have used permutation symmetry Kaf3 = Ki2. 
{N and are defined in Sec. ||.) The middle term of 



(4.30) can be simplified, since it is linear in TYmf which 
has the same value for any of the p-sublattice singlets 
that we are considering. The final result is 



'Hea{{Sa}) 



where 



K = 



^kY,{s^-Spf + c, 



and 



C 



2JsNS{2S - 1)2 



(4.32) 



(4.33) 



(4.34) 



By recalling that Js ^ N ^ (see eq. (3.5)) and that ^ 
S — Ns/p, it may be checked that ( 4.32 ) scales with N. 

In the limit N —^ oo, the operator becomes per- 
fectly defined and can be equated to a c-number, Snia, 
where riia is a unit vector denoting the classical direction 
of sublattice a. Then our effective Hamiltonian ( 4.32 ) 
simplifies further to 



(4.35) 

Earlier, Larson and Henleyi also obtained a second- 
order effective Hamiltonian starting from a different sort 
of mean- field theory: 





(4.36) 



Here rhi (called "z^" in Ref. |^) is the classical spin unit 
vector at site i and — \ Jijsrhj \ is the magnitude of 
the local field at i. Their result is valid for Hamiltonians 
more general than the "p-sublattice" category, on the 
other hand it is valid only in infinite systems. In the 
special case that sites are grouped into p sublattices, as 
assumed in this paper, it can be che cked that ( 4.36| ) does 
reduce to our effective Hamiltonian (4.35) for the N ^ oo 
limit, with Hq = (J12 — Jii)Ns. 



SEMICLASSICAL CALCULATION OF 
SINGLET ENERGIES 



The effective biquadratic Hamiltonian (4.32) is the ba- 
sis for an analytic approximation of the energy eigenval- 
ues for a p = 4 spin antiferromagnet, as detailed in Ref. ||. 



A further series of mappings is applied to the effec- 
tive Hamiltonian resulting in a semiclassical problem of 
one spin with a quartic Hamiltonian. Bohr-Sommerfeld 
quantization and tunneling considerations are then ap- 
plied to this one-spin problem to obtain energies and 
level splittings. The procedure here is identical to that 
in Ref. Q which worked out the triangular lattice case, 
except that the coefficients are different, being functions 
of Ji and J2 as well as the type and size of the lattice. 

This section recapitulates aiid somewhat extends the 
recent calculation by two of usEJ of the low-lying singlet 
energies of the four-sublattice Heisenberg antiferromag- 
nets, assuming long-range order in the N 00 limit. 
That derivation started from a four- spin system with the 
biquadratic effective Hamiltonian ( f.32 ) derived above. 
When the total spin is zero, as is the case here, the clas- 
sical dynamics is separable into a trivial and a nontrivial 
part (Sec. VA). The nontrivial part is very similar|_tp|_a 
single spin in an anisotropy field of cubic symmetry.t!Jt3 

Next, using Bohr-Sommerfeld quantization of the clas- 
sical orbits of this effective one-spin Hamiltonian, the 
low-lying singlet energies were shown to form clusters 
of two or three degenerate levels and the energies were 
approximately estimated (Sec. VB). It was found that 



the lower energy clusters are threefold degenerate, corre- 
sponding to excitations near the "coUinear" states of the 
classical ground state manifold (with four spins form- 
ing two pairs that point in opposite directions), and 
the higher energy clusters are twofold degenerate, corre- 
sponding to noncoUinear "tetrahedral" states (with four 
spins pointing at the four vertices of a tetrahedron) . 

Finally, tunneling between Bohr-Sommerfeld orbits in 
different energy wells generates splittings exponentially 
small (in the cluster size) within the level clusters, ex- 
cept when there is destructive interference of tunneling 
amplitudes associated with different paths which leaves 
the clusters unspiit. From the geometrical phase of the 
tunneling paths,E3 one could predict the pat tern of split 
and unspiit states in the clusters (Sec. V C ). The pre- 
dicted patteHi of level clustering was in accord with nu- 
merical datao for the J1-J2 triangular lattice with = 16 
and A^ = 28. 

We add a new discussion of the other quantum num- 
bers, derived from the spatial symmetries of the spin 
arrangements (Sec. VD). Since it is advantageous for 
exact diagonalizations to break the Hamiltonian into 
blocks with different symmetries, the symmetry eigenval- 
ues are often obtained as a byproduct, which may now 
be checked against our predict ions. 

The section concludes (Sec. VE) with a discussion of 
the system-size dependences of the various energy scales 
characterizing low-lying eigenstates. This illustrates how 
brute-force diagonalization of a (necessarily small) sys- 
tem can be misleading as to the energy scales of a large 
system, and how an analytic view can correct this. 
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A. Coherent states and separation of rotational 
variables 

We now carry out the semiclassical calculation on our 
quantum p-spin effective Hamiltonian (4.32) as in Ref. |. 
Given a quantum spin Hamiltonian, we can obtain up- 
per and lower bounds of the quantum partition function 
using the classical integral representation.cJ A classical 
spin Hamiltonian can then be derived_and is the start- 
ing point of semiclassical calculations.E3 We note there 
are multiple ways to define this semiclassical spin Hamil- 
tonian and two of which correspond to the upper and 
lower bounds of the partition function respectively (they 
are called Hq and Hp respectively in Ref. ^5|). Here we 
will use the Hq option which is obtained by taking expec- 
tations of the quantum Hamiltonian within spin coherent 
states. 

Using spin coherent states |{ma}) where riia is a unit 
vector (see, e.g., Ref. |6|), with 1 < a < p, we get, 



Heff({ma}) = ({rhallTioffKmc}) 



Ku ^ (m„ • m,, 

Q</3 



(5.1) 



where Ku = KS^{S - 1/2)2 ^ -pkS'^{AS{p - 

l)+p+l)/8 + C. For a closed orbit in this space {xa.a{t)} 
parametrized by the geometrical phase (important for 
semiclassical calculation) is, 



c& = ^5f]({m„(0}), 



(5.2) 



where f]({m(t)}) denotes the spherical area the trajec- 
tory m(t) has swept out on the unit sphere around the 
z axis. Because each of the p unit vectors lives on a 
two-dimensional spherical surface, so with the constraint 
^ riiQ, — 0, we have a {2p — 3)-dimensional space speci- 
fied by {ihc}. Furthermore, the global rotationjSeparates 
from the internal rotations (at least for p — 4;EJ this de- 
pends only on the permutation symmetry of eq. (4.32) ). 
We are left with a simpler problem on a two-dimensional 
spherical surface specified by a new unit vector n. 

From now on, we will work on the special case p — A 
which is the case for the fee and triangular lattices. First, 
let = (ihp -I- m4)/2, ^ = 1,2,3 (not a unit vector). 
It is easy to check that these three vectors are mutually 
orthogonal; therefore, given a coordinate system {e^}, 
we can rotate {n^}, with a proper rotation matrix TZ, 
so that 7?.n^ = n^e^, where are scalar coefficients. 
We then form a new vector fi = (ni, 712,713); it is easy 
to check that this is a unit vector (n^ -|- rij -|- rig = 1). 
In terms of this fi , w e obtain a one-spin semiclassical 
Hamiltonian from (5.1), 



Wi(n) 



3 



2Kr 



(5.3) 



It is easy to show that 72,mi = (n-i, — 712, — 713), 
7?.m2 = (-711,712,-713), TZm.3 = (-711, -712, rt3), TZrhi = 
(711,712,713), which says that each nia{t) traces out the 
same shaped trajectory as n{t) and makes the same con- 
tributation to the geometrical phase (5.2). Therefore, 
we have $ = 4S'r2(n(i)) = Sil{h{t)), where we have 
used S = AS = Ns. In the sense of semiclassical orbits 
and geometric al p hases, we have mapped the problem of 
Hamiltonian (|5.l| ) with four spins S to the problem of 
Hamiltonian ( |5.3| ) with one spin 5*. 

It must be pointed out that the definition of our co- 
herent states llrha}) included an averaging over all rota- 
tions. It follows that the mapping to n is not one-to-one. 
If TZjT is a TT-rotation about any coordinate axis , then 
(7?.7r7?.)n^ = n'pGp, and so the rotation matrix 7?.^ 7?. sat- 
isfies our condition just as well as TZ. In other words, in 
the n representation, two vectors related by a tt rotation 
about any coordinate axis represent the same four-spin 
state {rho,}. We need to take this discrete redundancy 
into account in the semiclassical calculation below. E3 



B. EBK quantization 

To obtain energies from the semiclassical spin Hamilto- 
nian, we use the EBK, or Bohr-Sommerfeld (BS), quanti- 
zation condition which says that the geometrical phase pf 
the Bohr-Sommcrfcld orbit is 2ttI where I = 0,1, 2SEB 
(The BS semiclassical calculation described in this sec- 
tion was carried out in Ref. |l] for the triangular system. 
The difference here for the fee system is the value of S*.) 
On the sphere, the spherical angle enclosed by two curves 
Oi{(f)) and ^2(0) is / (cos6'i — cos6'2)d0, where the integra- 
tion limits are determined by the geometrical locations 
of the two curves. For our quartic Hamiltonian, we first 
consider 



(5.4) 



where x = cos^ 9 and a(0) = cos^ (j) + sin'* </> = (3 
cos40)/4. Solving this quadratic equation, we get. 



cos^? = ±1 



, ± y^u ~ {1 — u)a 



(5.5) 



where the ± signs need to be chosen according to the 
location of the orbit. 

As mentioned in Ref. |l], two types of orbits on the n 
unit sphere are relevant. (1) "CoUinear" (C) states that 
orbit around the three coordinate axes (they are there- 
fore three- fold degenerate). The quantization condition 
for these states is $ = Anl (a factor of two is present be- 
cause of the rotation by tt redundancy mentioned above) . 
(2) "Tetrahedral" (T) states that orbit around two axes 
(±1,1,1)/a/3 (and equivalents) (they are then two-fold 
degenerate). Here the quantization condition is $ = 2^1. 

For the C-type orbits, the BS quantization condition 
leads to the equation, 
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'-i ^wn = 



(5.6) 



For T-type orbits, we get, 



— y/u— (1 — u)c 



da 



Vl-(4«-3)2 



where in the last integral we have changed integration 
variable to a, using d(t> = ±da/y/l - (4a - 3)2 (± signs 
determined by location). Without loss of generality, the 
C-orbit equation is obtained from orbits around (0,0, 1) 
and the T-orbit equation from those around (1, 1, 1)/V3. 
For given S and /, these two equations are solved nu- 
merically for u using Mathematica, and the semiclassical 
energy is then Kif{2 — 8u) + Cu. Semiclassical energies 
are presented in Tables III and |^ for fee and triangular 
lattices respectively. The result for the triangular case 
corrects a mistake in the last column of Table I in Ref . |l| 
which came from a wrong sign in the expression for Cu 
just under eq. (6) in Ref. The correct expression should 
be Cu - -52(55 + bl2)K + Cbiq. 



C. Tunnel splittings and eigenvalue fine structure 

The preceding section establishes that the low-lying 
levels are grouped into clusters of three in the lower part 
(corresponding to the three-fold degenerate C-type or- 
bits) and clusters of two in the higher part (correspond- 
ing to the two- fold degenerate T-type orbits). Average 
energies for the clusters have been calculated from the 
BS quantization condition. Finer level splittings within 
each cluster are results of tunneling between BS orbits, 
and degeneracies can remain due to destructive interffiCa 
ence of the geometrical phases of the tunneling paths. 
Using phase considerations, we can therefore obtain the 
split/unsplit pattern for energy clusters: it depends on 
N,s,l. 

We here simply quote the result from Ref. |l|. If we 
write "(i^i, 1^2)" for the pattern of eigenvalues in a cluster, 
meaning that the (lower, higher) eigenvalues have degen- 
eracies {vi,V2), respectively, then for the three- fold de- 
generate C-type orbits, we obtain a cluster pattern (2, 1) 
if the expression Ns/2—3l is odd or (1, 2) if it is even; and 
for the two-fold degenerate T-type orbits, we get (1, 1) 
when Ns j^Al is divisible by three, or (2), i.e., unsplit, 
otherwise.tj Notice that, ignoring whether the splittings 
are large or small, the pattern is 11221122... at either 
end. Indeed this sequence continues unbroken across the 
energy of the classical separatrix of the orbits, which is 
the boundary between the two clustering behaviors (i.e., 



where the C-type orbits meet the T-type) . As energy in- 
creases, the splitting within each cluster grows, until at 
the separatrix energy (it = 1/2) it is typically compara- 
ble to the separation between successive clusters, causing 
an ambiguity in labeling the levels. 

The case of fee type-I with N = 16 is trivial since 
SH = which means that the four-sublattice singlets 
are the exact ground states, and are exactly degenerate. 
Our calculation, of course, gets this right since it gives 
a zero coefficient in the effective Hamiltonian {K12 = 
4(iV - 4)jf/N and N = i)S5 Here, 5 = 2 and we have 
7)25 -f 1 = 5 degenerate singlet states (see Appendix ^). 

Table III shows the results For fee Type-I with N = 
32, s = 1/2. The preceding semiclassical considerations 
lead to an energy pattern which begins with (1,2)(2, 1) 
and ends with (1,1)(2). These clusters correspond to 
orbits Co, Ci, Ti, and Tq. The total pattern could be 
written as (1, 2)(2, 1/1, 1)(2), where the "1/1" indicates 
a single level which can equally well be assigned as the 
highest level in cluster Ci or the lower level of cluster Ti . 
In Table II], the entries with correspond to this level 



with the ambiguous labeling. The diagonalization results 
(especially for the full 32-site system, see Sec. VI C| below) 
illustrate how, near the separatrix energy, the inter- and 
intra-cluster splittings become comparable. 

Note that in this paper we have not attempted to 
estimate the tunnel splitting (i.e., the splitting within 
each cluster) quantitatively, because we have only com- 
puted the "WKB" exponent in this splitting and not its 
prpfartny fall thcsc Splittings are exponential in cluster 
size)MBEl 



D. Spatial symmetries 



The Hamiltonian (LJ) is invariant under the space- 
group £, which acts on the lattice of spins. Consequently 
the eigenstates transform as representations of this sym- 
metry group; in particular, they all (can be chosen to) 
have definite crystal momentum under lattice transla- 
tions. Similarly, our p-spin Hamiltonian (4.32) is invari- 
ant under all permutations of the p spins, so its eigen- 
states transform as the representations of the permuta- 
tion group Sp. „ 

Now, as observed by Lecheminant et alp each space- 
group operation induces a permutation of the p sublat- 
tices. Thus the symmetry eigenvalues of the p-spin eigen- 
states under the permutation group are the same as those 
of the corresponding eigenstates of the real Hamiltonian 
under the corresponding space operation. Ref. || com- 
bined this idea with group theory to count the number 
of times each representation appears among the p-spin 
singlets. Here, we outline how the same notion may be 
combined with semiclassics to identify which level has 
which quantum number, including the specific correspon- 
dence for both the triangular and fee realizations of p = 4. 
That will require an additional mapping from p-spin per- 



10 



mutations to (certa in) cubic-symmetry point operations 
on the n- vector in (|5.3|) . 

Our p-spin approximation groups the N spins of one 
sublattice into one big spin of maximum length, i.e. our 
states are totally symmetric under any permutation of 
sites within a sublattice. Thus we implicitly assumed 
that the symmetry eigenvalue is unity under any transla- 
tion that takes spins into the same sublattice. In effect, 
such translations are equivalenced to the identity, and 
the entire space group reduces to Sp (or a subgroup). 

In the four-sublattice case, every lattice translation is 
equivalent (if not to the identity) to a permutation of the 
class (12) (34). (We write permutations in the standard 
notation where each parenthesis is a cyclic permutation 
and each number is the label of a sublattice.) But for any 
configuration of the four classical directions {liia}, that 
permutation can be implemented instead by a 7r-rotation 
about the axis 113. Since we defined our coherent state 
llriiQ,}) to be symmetrized with respect to all rotations, it 
follows that the stat&_must be even under permutations 
of the (12)(34) class.E3 One corollary is that, in the 4- 
spin case, the low-lying singlets all have wavevector zero, 
as indeed rwas found in the exact diagonalizations of the 
triangulam and of the fee type 1 (Sec. VIC). 

Another corollary is that all permutations of this class 
are equivalenced to the identity; in effect, our permu- 
tation group is replaced hy S4/Q = S3 where Q is the 
subgroup consisting of (12)(34), (13)(24), (14)(23) and 
the identity. 

The permutations of class (123) correspond to the 
threefold axes in the triangular lattice or the fee lattice, 
as well as about the threefold axes (±1,1,1)/V3 of the 
n-sphere under Hamiltonian (5.3). Group theory suffices 
to identify the behavior of each state under (123): the 
nondegenerate states have eigenvalue 1 and the twofold 
degenerate states have eigenvalues e^^'^^/^. 

However, group theory is insufficient to identify the 
eigenvalues for the nondegenerate states under the odd 
permutations, which might be +1 or —1. The permuta- 
tions of class (12) correspond to any lattice reflection in 
the triangular lattice, or reflections in a {110}-type plane 
for the fee lattice. The permutations of class (1234) cor- 
respond to glide planes in either lattice, but they are 
equivalent to the (12) class modulo subgroup Q. On 
the n-sphere, an odd permutation is best represented by 
a 7r/2-rotations about one of the three coordinate axes. 
(Recall that the 7r-rotations are equivalenced to the iden- 
tity.) 

Consider the eigenstates in a three-cluster, mixed from 
three BS orbits Cix, Ciy, and Ciz (these are collinear 
states rotating around the X, Y, and Z axes respec- 
tively). The wavefunction described by Cix picks up 
a phase Airl when it goes full circlejJ thus it gets mul- 



rotation, or equivalently under the permutation (23). 
Thus, the lone eigenstate in the lowest cluster is even 
under odd permutations of sublattices, and alternates 
odd/even from then on. The twofold degenerate eigen- 
state contains can be broken into one odd and one even 
state (the state with eigenvalue —(—1)' under (23) is a 
mixture only of Ciy and Ciz)- 

Numerical diagonalization of the four-spin system (for 
all S < 13/2) co nfirm s the predicted pattern of even/odd 
states (see Sec. VIA ). In the exact-diagonalization re- 
sults, it can only be checked in the_case of the N — 
16 triangular lattice (which agrees). Ej Empirically, the 
even/odd alternation (of the nondegenerate states) con- 
tinues across the separatrix to the higher-energy clusters 
composed of two states, which are mixed from the two BS 
orbits called Ti±. It turns out that whenever such a pair 
is split, the lower level is even (odd) and the upper level 
is odd (even) under the odd permutations, according to 
whether I is even (odd). 



E. Dependence on size A'^ 

How do our results relate to the -macroscopic limit 
N ool It has long been understoocO that in quantum 
antiferromagnets symmetry-broken states (those with a 
nonzero expectation of each spin) are superpositions of a 
family of eigenstates which become degenerate with the 
ground state only as ^ 00. This motivates a classi- 
fication of the smallest energy gaps, measured from the 
ground state. 

In the present case of a four-sublattice system, the 
smallest gap is within a cluster of three singlets (as 
derived in Ref. |^) which are mixed from the lowest 
Bohr - Som merf eld orbits {Cxo, Cyo, Czo) as discussed in 
As discussed above, this gap is a tunnel 
I and is exponentially small in N. (The same is 



tipied by 



(—1) under a tt/2 rotation about the 



X axis (that orbits center). The nondegenerate eigen- 
state has an equal admixture of all three states, in par- 
ticular Cix, so it must have eigenvalue (—1)' under this 



Sec. |V] 
splitting 

true for any discrete symmetry breaking. But this gap is 
larger than in many other cases, since the tunnel barrier 
comes entirely from quantum fluctuations.) 

The other gaps are associated in a well-known fashion 
with the continuous symmetries or degeneracies of the 
classical system. If a continuous symmetry is broken and 
the interactions are short-range, there are always gapless 
"Goldstone mode" excitations, in this case the acoustic 
magnons. The states with one long- wavelength magnon 
(at a nonzero wavevector k) have much larger energies, 
E = 0(1 /N^/'^) as — > 00 in a d-dimensional system. 
This follows from the antiferromagnetic magnon disper- 
sion Lo = 0(|k|) and from |k| = 0{l/N^/'^) for the small- 
est wavevector (with periodic boundary conditions) . Our 
present theory does not address the magnon states, and 
their gap is not so small anyhow in a finite system. 

The limit of the acoustic magnon as k = has zero 
energy, but only in the N 00 limit. It is better to view 
it from the classical formula E = M'^/2Nx where M is 
the total magnetization and x is the susceptibility (per 
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spin) in the macroscopic hmit. Then provided Stot/N ^ 
1, the quantum state of total spin Stot has energy E = 
S'totl'S'tot + 1)/2A^X- This "tower" of multiplets is generic 
whenever the system breaks sp in , r o ta tion symmetry and 
the ground state is a singlet These states have 

total wavevector 0, but S'tot > 0. Their gap is much 
smaller than the magnons'; it is the next-smallest kind 
of gap after the tunnel splitting, lowest state of spin 1. 

Finally, just as the rotational symmetry breaking im- 
plies the Goldstone modes, the additional continuous de- 
generacy of our four-sublattice classical ground states im- 
plies a gapless "degeneracy mode" at harmonic order, 
with Wdcgcn(k = 0) = 0. However, quantum fluctua- 
tions - appro ximated by our biquadratic effective Harrail-. 
tonian ( 4.32 ) - break that degeneracy, creating a gapCJ'EJ 
'^dcgcn(k = 0) = A this gap becomes constant as A'^ ^ oo. 

On the other hand, A is exactly the energy spacing be- 
tween successive Bohr-Sommerfeld orbits surrounding a 
single collinear state on the n-sphere (orbits called "Cix" 
m Ref. 0). It can be verified that our semiclassical pre- 
scription indeed implies a constant limiting A in the large 
N li mit. Expanding the "one-spin" effective Hamiltonian 
(5.3) about a minimum point (e.g. the north pole), the 
l-th. orbit has energy Ui — Uo + 16Kij9f, where 9i is the 
radius of that (nearly circular) orbit on the unit sphere 
(assuming S = Ns is large enough that 6'; <C 1). Further- 
more the orbit ar ea i s fli — TrOf, so from the BS quanti- 
zation condition (5_^) we obtain [/; = C/q + 'A where the 
gap A = 6A{Ku/S). 

At this point we can trace back the N dependence 
of the parameters in the successive versions into which 
the Hamiltonian has been mapped (see Table ||). Writ- 
ing "J" to schematically represent all couplings, and 
"z" to represent the total coordination number, we get 
Js = OiJz/N), Ki2 = Oi.Pz), k = 0(J/8^3), and 
Ku = 0{JS/8) = 0{JS/32). Thus finally, A 
constant to leading order in 

The importance of the above discussion is that at 
N — 16, 28, or 32, the apparent A is still strongly size- 
dependent, and is smaller than the "tower" gap to the 
first eigenstate with spin M = 1. Only the analytic 
analysis reveals that the latter energy scale will become 
smaller in a large enough system. 



VI. NUMERICAL RESULTS FOR FCC TYPE I 

In this section, we present our numerical results from 
exact diagonalization for a lattice with various levels of 
appro ximation based on our effective p-spin Hamiltonian 
(4.32). First, we diagonalize it directly for p = 4 with 
parameters suitable for the fee Type I system. Then we 
will compare the predictions from four-spin diagonaliza- 
tion and semiclassics with an exact diagonalization of a 
32-site system. (The exact diagonalization of a 16-sitc fee 
Type I system gives five degenerate singl et st ates which 
are what we expect from theory (see Sec. VC).) Finally, 



0(2J), 



we will also give numerical results for A' = 16 and N = 28 
triangular lattices, correcting some minor numerical mis- 
takes in Ref. |l|. 



A. Diagonalization of a four-spin system 

We have exactly diagonalized the four-spin Hamilto- 
nian (4.32), as in Ref. |^, representing a fee lattice with 
iV = 32, s = 1/2 (TV = 8, 5 = 4). As discussed in Ref. ||, 
the classical antiferromagnetic Type I ordering occurs 
when Ji > 0, J2 < 0. Here we have used Ji — 1, J2 = —1. 
Our code was previously described in Ref. |lO|. In the di- 
agonalization here, we have used the conservation of total 



magnetization in the z direction, S^^^, to reduce the size 
of the Hamiltonian matrix, and also the (1234) cyclic 
permutation symmetry among the fo ur sp ins. (See the 
discussion of odd permutations in Sec. |VI]| .) We identify 
the singlet states as those which appear in the sector with 
total z spin 0, but not in the sector with total z s pin 1. 

The result is presented in Table III (and in Table [V for 
the triangular case). It should be emphasized that this 
calculation may be done for systems much larger than 
A^ = 32. We could handle 5* = 30 which corresponds to 
A^ — 240 for a spin 1/2, p = 4 system. 



B. Comparison to a one-spin Hamiltonian 

As we mentioned above, there are multiple ways to 
get a semiclassical spin Hamiltonian from a given quan- 
tum spin Hamiltonian. Different prescriptions agree to 
the leading order in spin, but the higher order correc- 
tions are much harder to obtain. To estimate the er- 



ror of our s emic lassical Hamiltonians Ticff ({iha}) (5.1) 



and 7ii(ri) (5.3) with respect to the quantum four-spin 
Hamiltonian 7icff({Sa}) (4.32), we construct a one-spin 
quantum Hamiltonian 7ii(S) using spin coherent state 
|n) = \e,(j)), such that (6*, 0|?i;i(S)|6', (/>) = Hi(n). (In 
the notation of Ref. 7ii(S) is the quantum Hamil- 
tonian for which Hi{h) is the Hq semiclassical Hamil- 
tonian.) 7ii(S) and Ticsii^a}) then correspond to the 
same semiclassical Hamiltonian 7Yi(n), and by compar- 
ing their eigenvalues we can estimate the error introduced 
in the semiclassical calculation. 

Here we need the expectation value of the quantum 
operator + Sy + in coherent states. We use the 
following identity (see Appendix 



St 



5(5-i)(5-l)(5-^) 



sin'' 9 sin^ ( 



Isis- 



1' 
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54 



Sy 



St~'^S']. (6.1 



Here we have emphasized quantum operators with S"^, 
and in the last Une, denotes classical spin variables 
{Sfj, — Sn^, where is defined in (^) above). We 
therefore get a one-spin quantum Hamiltonian from the 
semiclassical Hi (n) (5.3), 



where 



and 



Kn = 



S{S-l/2){S-l){S~?,l2y 



Co = ^ ( S\S + lf 
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(6.2) 

(6.3) 

-Cu- 
(6.4) 



Note that here S* = 45* is the total spin length. 

We have diagonalized this one-spin Hamiltonian (6.2) 
with 5 = 16 for the iV = 32, s = 1/2 foe system. It 
has the same semiclassical approximation as our four-spin 
system, except that half the eigenstates of the one-spin 
system are disallowed in the four-spin system. We em- 
phasize that the single-spin system is not a valid model 
for the extended A^-spin system in any approximation. It 
is included only to illustrate the errors in the semiclas- 
sical approximations mentioned above, since exact diag- 
onalization gives different results for the one-spin and 
four-spin quantum cases. Tables ^ and ^ con tain ex- 
act diagonalization results for this Hamiltonian (|6.2|) . 



of all spins; and working in subspaces with particular 
eigenvalues of these symmetries. The spin-flip symme- 
try is only useful for Sz — and the reflections are only 
useful if the translational eigenvalues e^^'^ correspond to 
wavevectors k = or k —{ti, 0, 0). For Sz — Q & subspace 
dimension of 1213429 was obtained, while for Sz — 1, 
where the spin-flip symmetry was not used, the dimen- 
sion was 2259363. 



Fortunately, the discussion in Sec. VD, shows that the 



predicted low-lying singlets all have k = and hence are 
included in the subspaces we checked. Note also that 
reflections of the {100} type map each sublattice into 
itself, i.e. they correspond to the identity permutation 
in 54 , so all the low-lying states are predicted to be even 
under this symmetry, which is confirmed by our results. 

The lowest lying eigenstates and eiaeii energies were 
calculated using the Lanczos algorithmoH with a ran- 
dom seed. If an energy level is degenerate, a single run 
of the algorithm will only find one eigenvector with this 
energy. In order to determine the degeneracy of the low- 
est levels we ran the algorithm several times, each time 
with a seed that was orthogonal to the low-lying states 
found in previous runs, so that only the remaining state 
space was investigated. Blue to numerical instabilities in- 
herent in the algorithm,c3 the orthogonalization was also 
applied after each step of the Lanczos procedure. The 
program was tested on a 16-spin fee cluster for which 
eigenvalues and degeneracies can be found analytically, 
and gave correct results. Further details of our method 
are found in Ref. ^. 

The exact diagonalization results for the fee lattice are 
in Table III. In Table [V we also show the results for 
the triangular case which was diagonalized in Ref^, cor- 
recting a factor-of-two mistake in Table I of Ref. |iP(The 
Hamiltonian used in Ref. |^ is twice the Hamiltonian in 
Ref. |l| and here.) 



C. Exact diagonalization of the A'^ = 32 system 

The lowest-lying states of a cubic 32-spin cluster with 
fee structure and periodic boundary conditions were cal- 
culated by numerical diagonalization. This is the small- 
est cluster that can accomodate all three types of antifer- 
romagnetic order predicted by mean field theory (AF-I, 
AF-II and AF-III), depending on the values of Ji and 
J2. It is a consequence of the periodic boundary condi- 
tions that for a given site, the six different next-nearest 
neighbor sites become pairwise identical, such that a spin 
is only coupled to three different next-nearest neighbor 
spins, each with the strength 2J2- Our calculations show 
that thig|dr)es not affect the type of order adopted by the 
system.B'E^ 

The size of the diagonalization problem was reduced 
by employing the following symmetries of the Hamilto- 
nian: total magnetization ^2, translations, reflections in 
planes of the {100} type, and the simultaneous flipping 



VII. CONCLUSION 

In summary, we have studied a class of antiferromag- 
nets such that their classical ground states have a p- 
sublattice structure and contain nontrivial degeneracies 
due to frustration. We are interested in the low-lying 
energy clustering patterns and have used a p-sublattice 
approach following the structure of the classical ground 
states. A degenerate family of singlet ground states is in- 
troduced, and an effective Hamiltonian is derived to ac- 
count for the fluctuations within sublattices and is used 
to calculate the low- lying energies of the antiferromagnet. 
This effective Hamiltonian couples the sublattice spins 
and is of the biquadratic form. It is written explicitly for 
a finite-size lattice and is therefore useful for comparison 
with exact diagonalization results. We diagonalized a 32- 
site fee Type-I system and compared its low-lying singlets 
with the analytically obtained eigenvalue patterns (from 
the effective Hamiltonian by diagonalization and by semi- 



13 



classical calculations). Our main analytical results are 
summarized i n T ables || and ^ an d our numerical results 
are in Tables III (for fee) and |IV| (for triangular lattice). 



In Table [V, for the triangular lattice with iV = 16 and 
28, we have corrected some numerical mistakes in Table 
1 of our earlier paper Ref. |l|, resulting in a much better 
agreement between exact diagonalization and semiclas- 
sical results. As seen in Tables [11 and IV, our method 
gives good agreement between theory and numerics for 
cluster splittings, but with an overall shift for absolute 
energy values. 

It seems likely that our general approach could be ap- 
plied to catalog the energy and symmetry eigenvalues 
in exact diagonalizations of some other lattice antifer- 
romagnets. By this approach, we mean (1) represent- 
ing the low-energy classical states by several sublattice 
spins, with lengths proportional to the size N; (2) writ- 
ing an effective Hamiltonian for the sublattice spins; and 
(3) analyzing this semiclassically, with especial attention 
to geometical phases. A caution is that the results are 
likely to be interesting mainly when there are multiple 
tunneling paths, and when the tunnel barriers are small 
compared to other energies (due either to classical degen- 
eracies broken by quantum fluctuations, or else to weak 
anisotropics) . 

The motivation for a careful study of the clustering 
patterns of low- lying eigenstates, and thus for this pa- 
per, is that the pattern might be diagnostic of the ul- 
timate long-range order, even in systems so small that 
correlation lengths and order parameters are inconclu- 
sive. Size is the key limiting factor to the usefulness of 
exact diagonalization for quantum many-body systems, 
so any approach that partially overcomes this is of inter- 
est. The eigenvalue pattern was decisive evidence of long- 
range order for the case of the s = 1/2 nearest- neighbor 
triangular antiferromagnet.E^ 

Even richer clustering patterns occur in systems with 
classical ground state degeneracies and/or discrete sym- 
metry breakings, as treated in this paper. However, stud- 
iesli3 subsequent to Ref. |] cast this partially into doubt: 
the eigenvalue pattern of the four-spin problem with full 
rotational symmetry reappeared for three spins with XY 
symmetry, and presumably in many systems with a dis- 
crete threefold symmetry (see Sec. V A of Ref. 10|). 

An alternate application of this theory would be to the 
burgeoning topic of small magnetic molecules, the mag- 
netic Hamiltonian of which may often be approximated 
as a combination of several large subspins.E^ In the ex- 
amples studied to date, the intramolecular interactions 
were ferromagnetic and highly anisotropic, but it is quite 
plausible that a molecule (or complex) with high sym- 
metry could realize a Hamiltonian similar to one of those 
treated in this paper. If furthermore an experimental 
probe became available to measure a large number of the 
energy levels, this subject could be developed as chemists 
have developed the semiclassical treatment of moleculat 
rotational-vibrational levels over the past twenty years.E3 

It would be desirable to generalize our derivation of 



the effective Hamiltonian to other systems with classical 
degeneracies, such as the fee Type II or Type HI antifer- 
romagnets (for comparison to the exact diagonalizations 
of Ref. ^). This does not appear trivial. We started by 
finding the quantum ground-states of the infinite-range 
Hamiltonian TYmf, and the degeneracy of its ground 
states was essential to our derivation, (the p-sublattice 
singlets) particularly in factoring out the energy denom- 
inators in Subsec. IV C . p 

Consider, as the simplest exampleH, the square lat- 
tice with sufficiently strong antiferromagnetic J2 that the 
classical ordering wavevector is (tt, 0) or (0, tt). The even 
spins form two sublattices, oriented oppositely, so they 
combine to make a net singlet. Likewise the two sublat- 
tices of odd spins make a singlet, so T^mf has a unique 
ground state in this case. (The classical degeneracy is 
reflected only in the excited states of TYmf-) Rather than 
pursue such a complicated derivation, one could instead 
just derive the effective Hamiltonian as in Ref. I (see our 
Eq. ( 4.36| )), but this does not capture the system-size 
dependence of the coefficients. 
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APPENDIX A: ADDING SPINS TO GET 
SINGLETS 

To add up p spins of length S to get a singlet, we use 
the well-known decomposition of the product of two rep- 
resentations: V^^ X V^^ = pISi-SsI ^^pISi-Sal-Hl _^ _^ 
X)Si+S2 ^ where T>^ is the spin S irreducible representation 
of SU{2). (Note that to get a singlet (P") in the preceed- 
ing expression, we must have 5*1 = 6*2.) Here we show a 
somewhat nontrivial calculation for p = 5 (and assume S 
is an integer). One can either find the coefficent ofD^ in 
V'^ xV^ xV^ X V^, or try (V'^ xV^ x V^) x [V^ x V^). 
Adopting the latter approach, and denoting Ns the co- 
efficent of in T>^ xV^ xV^, we get Ns = 2S+1 for 
< S" < 5 and iVs = 35" - + 1 for ^ -f 1 < 5 < 3^. 
The number of singlets for p = 5 is thus X]s=o -^s ~ 
5S{S + l)/2 + 1. (The upper limit of the sum goes to 
2S because the maximum spin length from x T)^ , the 
last two terms in the product of five terms, is 25'.) 



APPENDIX B: QUARTIC OPERATOR AND 
SPIN COHERENT STATE 

Here we note briefly how we compute the spin coherent 
state expectation value in (6.1). The straightforward way 
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is to use a generating function (Eq. (6.22) in Ref . [26|) . We 
show here a second way using tensor operators (Ref. 
Eq. (8Q)): 



(e,0|3^LM(S)|0,<^) 



(25)! 



{2S-L)\2^ 



LM( 



(Bl) 



Here Ylm{Oj 0) is the usual spherical harmonic function. 
3^lm(S) is the spherical tensor operator obtained from 
the polynomial S^Ylm{0, </>) by first substituting the op- 
erators Sx, Sy, and for classical numbers S^, Sy, and 
Sz, and then symmetrizing the resulting operator, e.g., 
the polynomial SxSy becomes the symmetrized operator 

{SxSy + SxSy) /2. 

This identity tells us that if we can write a spin opera- 
tor as a linear combination of spherical tensor operators 
3^lm(S), then its coherent state expectation value is easy 
to obtain. For our problem, it is easy to check the fol- 
lowing classical identity: 



Sx + Sy + ■ 

^^4,o(^,</') + 

5 



3 C.4 _ /47r 4 

{Yi49,(b)+Yi^^i{e,(b)) 



(B2) 



Then according to (Bl), we have 




-3^4,o(S) 



54 



70 

s ^ 



(3^4,4(8) +3^4,-4(8)) 



St - ^S' 
o 



(B3) 



What is the operator on the left hand side of the previous 
equation? Using the table in Ref. |4^, it can be checked 
that this spin operator is nothing but 



S S^, S y 



S^{S + lf -\s{S + l) 



(B4) 



On the other hand, we see that the three spherical har- 
monic functions involved here all have L — -i and ac- 



cording to (Bl) the coefficient in front of the expecta- 
tion value depends on L only, i.e., the coherent state 
expectation values of the corresponding spherical tensors 
have the same coefficients. We can then start from the 
classical expectation value 5*^ 



■ Sy + S^ 



3 c4 



S and find 



the corresponding operator. It is straightforward to get 
Sx + Sy + 5*2, but because of operator symmetrization, 
54 goes to S^iS + 1)2 - S{S + l)/3 (see Ref. ||, p.268, 
also Ref. p. 31 ). P utting these two terms together, we 



get the identity (6J) in our paper 
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TABLE 1. Parameters for the effective Hamiltonian 



General 



Triangular 



Fee 



P 



zip (« / /?) 

zip {a / /?) 



■^12 = (E„ Jr^zr2)/N 

Emf = (p/2)(JnS2 - Ji2S{S + 1)) 



4 
2 
2 



2(Ji + J2)/iV 


-4(Ji + J2)S'(S + l)/iV 
2( Jf + J|) - 4( Ji + Jif/N 



4 
4 


6 

4Ji/N 
6J2/N 

-8JiS{S +1)/N + 12J2S^/N 
4{N -4)Ji/N 



TABLE IL The series of Hamiltonians 



No. of spin Spin 



Hamiltonian 



Parameters 



Eqn. No. 



N = N/p S^Ns n,ff{{S^}) ^ ~Kj2^^p{Sc ■ Spf + c 
1 S^Ns m (n) = -8Ku E^=i K + 2Ku + Cu 

1 s M^)--KQY.UJt + CQ 



J m 

Js ^ Ji2~N.hi/{N -1) 
K = Ku/{2JsNS{2S - if) 
C = -pkS^{{p + 1)S^ - 2{p - 1)S - 2)/2 
Ku = KS^{S~ l/2f 
Cu = -pKS'^iASip -l)+p+ l)/8 + C 
Kq = 8Ku/{S{S ~ 1/2){S - 1){S ~ 3/2)) 
Cq = 3 (5^(5 + if ~ Sis + l)/3) Kq/5 - UKu/5 + Cu 



(1.1) 



(4.32) 



(|5J 



TABLE IIL Eigenenergies AE ^ E - Emf for iV = 32 fee Type I ( Ji = 1, J2 = -1). (Entries denote a level where 

cluster labeling ambiguity occurs and are explained in Sec. [V C|.) 



Orbit, degeneracies 



Exact(l.l) 



(mean, levels) 



4-spin ( |4.32D 
(mean, levels) 



1-spin (|6J) 
(mean, levels) 



Semiclassic (5.J) 
(mean) 



To 2 

HI 

Ci 



-1.5388 
-1.7457- 

-1.8409 • 

-2.3178 ■ 



-1.7126 

-1.7788 

-1.8720 
-2.3151 

-2.3233 



-1.0316 
-1.2029 • 

-1.2849 • 

-1.6706' 



-1.1492 

-1.2566 

-1.2990 
-1.6703 

-1.6713 



-1.0659 
-1. 



-1.2519 
-1.5640 



-1.1436 

-1.2361 

-1.2599 
-1.5639 

-1.5642 



-1.0727 
-1.1718 

-1.2899 

-1.5639 



TABLE IV. Eigenenergies AE = E - Emf for TV = 16, iV = 28 triangular lattices (Ji = 1, J2 = 0.7). 



Orbit, degeneracies 


Exact(l.l)i 


4-spin (4.32) 


1-spin (3.2) 


Semiclassic (5.3) 




(mean, levels) 


(mean, levels) 


(mean, levels) 


(mean) 



TV 
To 2 

Co 



16 



{I 



TV = 28 
To 2 

2 



Ci 



-0.0756 
-0.1359. 




-0.1331 
-0.1417 



-2.0022 

-2.0275 
-2.3840 

-2.3845 



-0.0730 
-0.1366' 




-0.1343 
-0.1412 



-2.2993 

-2.4097 
-3.1539 

-3.1595 



-0.0792 
-0.1178' 




-0.1170 
-0.1193 



-2.2601 

-2.3261 
-2.9246 

-2.9262 



-0.0824 
-0.1176 

-2.0135 
-2.3699 

-2.9252 
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